from scipy.interpolate  import lagrange

def ployinterp_column(s, n, k=1):
    y = s.reindex(list(range(n - k, n)) + list(range(n +1, n+1+k)))
    y = y[y.notnull()]
    return int(abs(lagrange(y.index,list(y))(n)))